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Abstract 

We investigate the energy dependence of the astrophysical S factors for the 
reactions 7 Be(p, 7) 8 B, the primary source of high-energy solar neutrinos in 
the solar pp chain, and 16 0(p, 7) 17 F, an important reaction in the CNO cycle. 
Both of these reactions have predicted S factors which rise at low energies; we 
find the source of this behavior to be a pole in the S factor at a center-of-mass 
energy E = —Eb, the point where the energy of the emitted photon vanishes. 
The pole arises from a divergence of the radial integrals. 
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The 7 Be(p, 7) 8 B reaction, at center-of-mass energies E near 20 keV, plays an important 
role in the production of solar neutrinos 0. The neutrinos from the subsequent decay of 8 B 
provide the high energy neutrinos to which many solar neutrino detectors are sensitive. The 
cross section for this reaction is conventionally expressed in terms of the Syj factor, where 
the S factor is defined in terms of the cross section a by 

S(E) = a(E)Eexp(27ir](E)) , (1) 



where rj(E) = ZiZ^otJ yuc 2 /2E is the Sommerfeld parameter for nuclei of charges Z\, Z2 and 
reduced mass fi. The exponential factor in the definition of S removes the rapid energy 
dependence of the cross section due to Coulomb repulsion between the two nuclei. In the 
stellar core the probability of capture of protons by 7 Be, obtained by folding the thermal 
distribution of nuclei with the cross section, peaks at ~ 20 keV. Because the cross section 
diminishes exponentially at low energies, the only method of obtaining information about 
S17 at those energies is to extrapolate data taken at experimentally accessible energies (E > 
100 keV). 

The 16 0(jo,7) 17 F* (| + ; \, 0.498 MeV @) reaction, which occurs in the CNO cycle, is of 
little importance for energy production in the sun but of greater importance for hotter stars. 
As is the case for 7 Be(jo, 7) 8 B, extrapolation of data taken at high energies is necessary to 
obtain the S factor at energies applicable in the stellar core, E ~ 25 keV. 

Direct capture calculations of these two reactions predict an upturn in the S 
factor at threshold. As the capture in both reactions is primarily external, the S factors 
at astrophysical energies are determined by the product of the spectroscopic factor, the 
asymptotic normalization of the final (bound) state wave functions, and a purely Coulombic 
term. As the spectroscopic factor is independent of energy, the energy dependence of the 
S factor, away from resonances, may be studied without detailed knowledge of the nuclear 
structure. 

In each reaction the weakly bound final state (E-q = 137.5 keV [0] for 8 B and E-q — 
105.2 keV for the first excited state of 17 F 0) causes the S factor to rise as the center-of- 
mass energy approaches 0. In the case of 16 0(p, 7) 17 F* both the data and direct-capture 
calculations |4]|| show clear evidence of this low-energy rise. The 16 0(p, 7) 17 F capture to the 
ground state shows no such rise because the final state is more deeply bound and has higher 
angular momentum. For 7 Be(p, 7) 8 B the upturn occurs below the lowest experimental point 
so it is only observed in the calculations. 

Williams and Koonin || do an explicit expansion about zero energy for S 17 and give the 
first two coefficients in a Taylor series for the logarithmic derivative of Sn as —2.350 MeV -1 
and 28.3 MeV~ 2 . Using these in a (1,1) Pade approximant gives: 

(1 + 4.85E) 

" (1 + 7.20E) { ) 

where E is in MeV. This Pade approximant has a pole at —139 keV which is very close to 
their bound state energy of 136 keV. Thus we see that in the region of the threshold the bound 
state is important and induces a pole. In fact, a Taylor series expansion would converge 
only with a radius of the binding energy - - barely to the region that is experimentally 
accessible. Hence any functional form for the extrapolation to zero energy should contain 
the contribution from the pole. 
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Following ||, we write the astrophysical S factor as 

S = C(I 2 + 2I 2 2 )E* (j u ^ + J 12 ft 2 ) l _ , (3) 

where 

h = dr r 2 ifj iL (r)if) f (r)/k (4) 

J 

In Eq. (0), is the spectroscopic factor for a given angular momentum, L, and channel 
spin, S, (3ls is the asymptotic normalization of the bound state wave function, E~ is the 
photon energy, and k is the momentum of the incident proton. The final bound state wave 
function ipf(r) is normalized asymptotically to ipf{r) = W a j(Kr)/r while the initial wave 
function reduces to the regular Coulomb wave function divided by \J2nrj/ (e 27Tr) — 1). The 
unusual choice of normalizations is just to eliminate uninteresting factors from quantities of 
interest. Most of those factors have been collected in the coefficient C. 

To investigate the behavior of the integrals in Eq. (^), we first consider if)f(r) = 
W a ^(nr)/r for all radii and take ipio(r) = F (kr) / {^2nri/ (e 27T ' n — 1)). The integral for 
the s-wave then becomes 

roc i 

I = drrW a)1 (Kr)F (kr)/(k^2Trr]/(e 27TV -l)). (6) 

J 

At threshold, the integrands are peaked at large r: 40 fm for 7 Be(p, 7) 8 B, and 65 fm for 
16 0(p, 7) 17 F*. The tails of both integrands extend well beyond 100 fm, and are, in each case, 
indicative of halo states. The integral is smooth as k passes through zero and diverges as 
k — > in (E — > —Eb)- The nature of the divergence is determined by the asymptotic forms of 
the Coulomb wave function and Whittaker function for large r. For large r the Whittaker 
function is proportional to r -\ r i k \/ K e - Kr || (jjk is independent of k). While above threshold 
the Coulomb wave function oscillates at large radii, below threshold it is exponentially 
growing and is proportional to r' r? 'e' /c ' r . Thus the integrand approaches 

r .i-W(i/»-i/l*l) e xp[-(«- |jfe|)r] (7) 

for large r and the integral diverges as 

h oc 1/{k - \k\) 2 oc 1/(E B + E) 2 = 1/E 2 . (8) 

Since the integrand diverges as 1/E 2 the leading term and first correction term are both 
determined purely by the asymptotic behavior of the wave functions. The first correction 
term is not simply 1/E 1 but also involves logarithmic terms coming from the r~' r?fc !( 1 / K ~ 1 /' fc !) 
factor. The second correction term, of order E®, is not determined purely by the asymptotic 
value of wave function alone but also depends on the wave function at finite r. 

From Eq. (||), we see that the quadratic divergence of I gives rise to a simple pole in S 
at E 1 = 0. This suggests writing the S factor as a Laurent series: 

S = d^E^ 1 + d + d x E^ + ... (9) 
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As before the coefficients of the first two terms, d-\ and do, are determined purely by the 
asymptotic forms of the wave functions while the third coefficient, d%, is also dependent on 
the short range properties of the wave functions. 

In Fig. m, we present the data of Morlock et al. || for the 16 0(p, 7) 17 F* S factor (top) 
and for the product E 1 S (bottom). In the top panel, the energy dependence of the S factor 
is well approximated by the form: 

^ = n — Ei = n ~u< ET" (10) 

Ey E + E-q 

where the constants c% and n are determined by the straight line fit to E^S shown in the 
bottom panel. The numerical values are given in Table |. There is remarkable agreement 
with the data except near the resonance at 2.504 MeV. Eq. (|T0"D is a convenient form for 
fitting experimental data and is motivated by both the Pade approximant and Eq. ([|). 
In Fig. |2|, the data of Filippone || (circles) and Kavanagh [ TDj] (diamonds) for the S 



factor (top panels) and the product E y S (bottom panels) for the 7 Be(p, 7) 8 B reaction are 
presented for energies well below the E = 633 keV Ml resonance. We take the data as 
normalized by Johnson et al. JTT[] to = 157 mb. The curves are similar in form to 
Eq. fllP]) , but with a quadratic term added, 



1 + Cl E + c 2 E 2 l + Cl E + c 2 E 2 
b = n = n . (11) 

Ey E + E-q 

The values of n and q for this reaction are also listed in Table |[ Different values of the 
normalization n are required to reproduce the Filippone and Kavanagh data, but the c, are 
determined from the threshold energy dependence of a direct-capture calculation following 
Ref. IJ. A cut-off radius of r = 2.3 fm was chosen to be consistent with the phase shift 
and energy dependence found by Barker ||12||. The upturn at threshold is clearly observed 



in the results of the calculation. The data are insufficient to determine this behavior or, 
equivalently, q. It will be very difficult to experimentally confirm this upturn since it is only 
pronounced below 100 keV. Fortunately it is theoretically well understood and both n and 
Ci depend primarily on the asymptotic normalization, spectroscopic factor and properties of 
the Coulomb force. Note that the curves presented in Fig. |2] should not be mistaken for a 
serious attempt at determining the Syj factor at zero energy; rather, they are illustrative of 
the energy dependence. 

The straight line approximation for E 7 S is valid for the 16 0(p, 7) 17 F* reaction up to 
~ 3 MeV. However, the quadratic approximation for the 7 Be(p, 7) 8 B reaction is not valid for 
energies above 0.4 MeV. Initially, the breakdown is caused by the resonance at 0.633 MeV. 
Above the resonance higher order terms in E, arising predominantly from <i-wave direct 
capture, become significant. 

In conclusion we see that the threshold peak in the S factor is associated with weakly 
bound states and arises from a pole at E^ = 0. Those bound states in 8 B and 17 F are halo 
in nature and so the associated radial integrals are, by necessity, long range. 

The authors would like to thank F.C. Barker, L. Buchmann, H. Fearing, T. Hemmert, 
K. Langanke and E. Vogt for useful discussions and R. Morlock for making tables of his data 
available. Financial assistance from the Natural Sciences and Engineering Research Council 
of Canada is gratefully acknowledged. 
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FIG. 1. Astrophysical S factor (top) and E 7 S (bottom) for 16 0(p, 7) 17 F*. The data of Morlock 
et al. [§] are compared to the fit as described in the text (solid line). 
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FIG. 2. The low-energy part of the astrophysical S factor and E~S for 7 Be(p, 7) 8 B. The data of 
Filippone |J (circles) and Kavanagh (n| (diamonds) are compared to the results of the calculations 
as described in the text (solid line). 
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TABLES 

TABLE I. The numerical constants q and E-q used to determine the energy dependences, with 
the normalizations n used when displaying the curves with the discussed data sets. 



Reaction 


ci (MeV" 1 ) 


c 2 (MeV~ 2 ) 


E B (MeV) 


n (keV 2 b) 


16 0(p, 7 ) 17 F* 


1.18 





0.1052 


1.59 x 10 3 


7 Be(p,7) 8 B 


5.36 


1.80 


0.1375 


3.74 a 










2.99 b 



(a) Kavanagh et al. [|10| 

(b) Filippone et al. |J 
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